Density Induced Interchange of Anisotropy Axes at Half-Filled High Landau Levels 



(N 
O 

o 

(N 



c3 



c3 



i 

C 

o 
o 



> 

o 

m 
O 
(N 
O 

-i— > 
03 



i 

C 

O 

o 



J. Zhu, 1 W. Pan, 2 H.L. Stormer, 1 ' 3 ' 4 L.N. Pfeiffer, 4 and K.W. West 4 

'Department of Physics, Columbia University, New York, New York 10027 
2 Department of Electrical Engineering, Princeton University, Princeton, New Jersey 08544 
3 Department of Applied Physics and Applied Mathematics, 
Columbia University, New York, New York 10021 
* Bell Labs, Lucent Technologies, Murray Hill, New Jersey 07974 
(Dated: February 1, 2008) 

We observe density induced 90° rotations of the anisotropy axes in transport measurements at 
half-filled high Landau levels in the two dimensional electron system, where stripe states are proposed 

=9/2, 11/2, etc). Using a field effect transistor, we find the transition density to be 2.9 x 10 n cm -2 
at i/=9/2. Hysteresis is observed in the vicinity of the transition. We construct a phase boundary 
in the filling factor-magnetic field plane in the regime 4.4 < v < 4.6. An in-plane magnetic field 
applied along either anisotropy axis always stabilizes the low density orientation of the stripes. 

PACS numbers: 73.40.Hm 



The physics of electrons in higher Landau levels has re- 
cently drawn much attention due to the discovery of large 
transport anisotropies at half fillings (v = 9/2, 11/2 etc) 
in high quality two dimensional electron gases (2DEG) 
in GaAs[]l], g]. At these filling factors, longitudinal mag- 
neto resistances R xx of samples grown on (OOl)-oriented 
GaAs show strong maxima in the (110) direction and 
deep minima in the (110) direction. The ratio of the re- 
sistance values can exceed 1000:1 in high quality samples. 
This large anisotropy is viewed as strong evidence for the 
formation of a unidirectional Charge Density Wave or so 
called "stripe state" around these filling factors || ||, ||, ||. 
Calculations based on the Hartree-Fock approximation 
show the stripes to form as the result of the competi- 
tion between the short-range attractive exchange inter- 
action and the long-range Coulomb repulsion. Recently 
there appeared calculations beyond the Hartree-Fock 
approximation 0, [§[ ^, [l0| |xi[| . Some of them propose 
the existence of liquid crystalline states with stripe order- 
ing and broken rotational symmetry 0, [hJ. Employing 
an edge state transport mechanism, the "easy" (low resis- 
tance) axis is along the stripes whereas the "hard" (high 
resistance) axis is perpendicular to the stripes. This iden- 
tification is supported by experiments JT^] . 

In spite of the strong evidence supporting the forma- 
tion of a stripe state, the origin of its preferred orien- 
tation remains poorly understood [fl3|, ft4| , 15 1. Theory 
does not a priori provide a preferred direction for the 
stripes. However, all previous experiments have iden- 
tified the (110) direction as the "easy" axis. It is often 
assumed that anisotropic imperfections at GaAs / AlGaAs 
interface act as the native symmetry breaking potential. 
Anisotropic roughness has been observed in GaAs layers 
grown by Molecular Beam Epitaxy [|l6| . Two groups have 
investigated the correlation between the surface rough- 
ness and the orientation of the stripes which form a few 
thousand A below the surface. Atomic Force Microscope 
images reveal the presence of roughness elongated along 



both the (110) and the (110) directions. However, its 
correlation with the orientation of the stripes remains 
controversial |lj| |l8). Theoretical studies on the influ- 
ence of a simple periodic modulation on the orientation 
of the stripes suggest that they prefer to align perpen- 
dicular to a weak potential JToj, p0[ . However, a parallel 
alignment may occur in the strong potential limit p(i|. 

Previous experiments that identified the (110) direc- 
tion as the "easy" axis were limited to fixed densities 
in the range of 1.5 - 3.0 x 10 u cm _2 §], § @. Fur- 
thermore, the correlation between the orientation of the 
stripes and the native symmetry breaking potential is 
difficult to deduce since the latter is likely to vary from 
specimen to specimen. Employing a tunable density 
Heterostructure Insulated Gate Field Effect Transistor 
(HIGFET) in our experiments, we have the unique op- 
portunity to access a wide density regime and to change 
the 2DEG density in situ without affecting other param- 
eters. Our data show that as the density of the 2DEG is 
raised above 2.9 x 10 1:L cm _2 , the "easy" axis rotates from 
the (110) direction to the (110) direction. This result 
demonstrates that the pinning direction of the stripes is 
not unique. 

Our HIGFET consists of a 600/im square mesa with 
edges along the (110) and the (110) directions. The 
structure of the specimen consists of a (001) GaAs sub- 
strate, overgrown by MBE with 0.5^m of GaAs. A 50A 
layer of AlAs is deposited, followed by a 4000A layer of 
Alo.32Gao.68As and capped by a heavily doped GaAs n+ 
layer, serving as a top gate. The AlAs layer at the inter- 
face with GaAs is intended to improve mobility. Sixteen 
Ni-Ge-Au contact pads are placed evenly along the edges 
of the mesa using standard optical lithography. One cor- 
ner pad provides the contact to the top gate (see in- 
set to Fig.|l|). Biasing the gate, we are able to contin- 
uously change the 2DEG density from 5 x 10 9 cm~ 2 to 
4.9 x 10 n cm~ 2 . The peak mobility is 1.1 x 10 7 cm 2 /Vsec 
at 2.3 x 10 11 cm -2 . The 2DEG density has a reproducible 
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FIG. 1: Interchange of anisotropy axes of the v — 9/2 state 
with increasing densities. The two central insets show the 
contact configurations used to measure the anisotropy. Data 
from the (110) configuration are represented by solid lines and 
data from the (110) configuration by dash-dotted lines. In the 
left panel, the stripes align along the (110) direction. In the 
right panel, they align along the (110) direction. See insets. 
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FIG. 2: Hysteresis at 2.6T in the (110) contact configu- 
ration^) and the (110) contact configuration(b). Sweep di- 
rections are indicated by arrows and the hysteretic areas are 
hatched. Sweep rate is ~ 1.0 x 10 10 cm~ 2 /min. c) Area of 
the hysteresis loop as a function of magnetic field, d) Filling 
factor range of the hysteresis in the ^-B plane. Symbols cor- 
respond to (a). The solid line represents the general position 
of the phase boundary. 



linear dependence on gate voltage. Only the lowest elec- 
tronic subband is occupied for the whole density range. 
Fragile fractional quantum Hall states such as the 5/2 
state and those around 3/2 are well developed for a wide 
range of densities, attesting to the high quality of the 
sample. The anisotropies at 9/2 and 11/2 are large, en- 
abling us to identify the anisotropy axes unambiguously. 

The electrical measurements are carried out in a dilu- 
tion refrigerator with a base temperature of 20mK using 
standard low frequency (13-21Hz) lock-in techniques. An 
excitation current of lOnA is used to avoid electron heat- 
ing. Ra;a; traces are taken at fixed perpendicular mag- 
netic field while sweeping gate voltage to change density. 
The left and right panels of Fig. |l| show traces taken at 
2.1 Tesla and 3.1 Tesla, respectively. At 2. IT, the 9/2 
state occurs at 2.3 x 10 n cm~ 2 . At 3. IT, it occurs at 
3.4 x 10 n cm~ 2 . The two central insets in Fig. [jjshow the 
two orthogonal contact configurations used to measure 
the anisotropy pH. They are labelled (110) and (110) ac- 
cording to the direction of the current flow. The (110) 
configuration is always represented by a solid line and 
the (110) configuration by a dash dotted line. In the 
left panel, at 2. IT, the "easy" axis is along the (110) di- 
rection, in agreement with all previous experiments. In 
the right panel, at 3. IT, the "easy" axis is rotated 90° 
and points now along the (110) direction. We conclude 
that the underlying stripe state has also rotated 90° and 
is now aligned along the (110) direction. This demon- 
strates that the pinning direction of the stripes is not 
unique. This observation is the most striking result of 
our experiments. 



Our results indicate that the pinning mechanism of the 
stripes is more complex than previously assumed. A few 
conclusions can be drawn. First, contrary to earlier re- 
sults, a (110) orientation of the stripes does exist in a 
perpendicular magnetic field. Second, assuming that the 
stripes are pinned by a unidirectional periodic potential, 
these stripes can be aligned either parallel or perpendic- 
ular to it. Although we do not know the mechanism 
responsible for the reorientation, we identify several con- 
sequences of changing density and discuss their effects. 
With increasing density, the electron wave function is 
squeezed and pressed harder against the interface. There- 
fore, the electrons experience the interface potential more 
strongly. On the other hand, with increasing density, 
screening also increases, effectively weakening the inter- 
face potential. In addition, the period of the stripes, 
calculated to be a few magnetic lengths at v = 9/2 [51 , 
decreases with increasing density. Any of these factors 
could lead to the observed reorientation of the stripes [p0|. 
Furthermore, we can not rule out the possibility of the 
coexistence of two orthogonal periodic potentials in our 
sample. Changing the period of the stripes could shift 
the relative effectiveness of these two potentials and may 
cause a reorientation. 

To find the transition density for the reorientation of 
the stripes, we take density sweeps at fixed magnetic 
fields from 2.0T to 3.3T in steps of 0.1T and identify 
2.7T and the corresponding density n= 2.9 x 10 11 cm -2 
as the transition point at v — 9/2. From 2.2T to 3.0T, we 
observe hysteresis between density up-sweeps and down- 
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sweeps in a narrow range of filling factors (Av ~ 0.2) 
around 9/2. As an example, Fig. ||a and b show the 
hysteresis at 2.6T and 40mK in both contact configu- 
rations. Sweep directions are indicated by arrows and 
the hysteretic area is hatched. The hysteresis decreases 
with increasing temperature and disappears altogether 
at 90mK while the anisotropy remains although it is 
weaker. These 90mK data unambiguously identify the 
orientation of the stripes in close vicinity to v = 9/2. 
Fig. ||c shows the area of the hysteresis loop in the (110) 
configuration at 40mK as a function of magnetic field. 
It peaks near 2.5T and vanishes at 2. IT and 3.0T. In 
Fig. ||d, we plot the range of the hysteretic behavior in 
the i^-B plane. This hysteresis "ellipse" is bounded by 
v = 4.4 and v = 4.6, which coincides with the filling fac- 
tor range in which the stripe state is expected to exist 
theoretically 

In general, the presence of hysteresis indicates a first 
order phase transition. As the system crosses the tran- 
sition point, the initial phase continues to survive due 
to some pinning mechanism. In our sample, the phase 
transition corresponds to the reorientation of the stripes 
while impurities act as the pinning mechanism. In the 
2/-B plane, the stripe state exists in a horizonal band 
between v—\A and v=A.Q. The hysteresis "ellipse" is 
embedded within this band, indicating the transition re- 
gion. To the left of the "ellipse" , the stripes have the 
(110) orientation. To the right of the "ellipse" , the stripes 
have the (110) orientation. Therefore, the boundary sep- 
arating the two phases must stretch through the whole 
bubble from the left side to the right side. We know 
from the 90mK data that the boundary intercepts the 
v =4.5 line at 2.7T. Furthermore, the orientation of the 
stripes along the border of the hysteresis "ellipse" can 
be deduced from data such as Fig. ||a and b since hys- 
teresis always prolongs the initial phase. For example, 
at 2.6T, on density up-sweeps, the stripes are along the 
(llO) direction, as seen from the minimum in the (110) 
configuration and the maximum in the (110) configura- 
tion. This indicates that the stripes near v = 4.4 have 
the (110) orientation. Conversely, density down-sweeps 
indicate the stripes near v = 4.6 have the (110) orienta- 
tion. Similar identifications can be made for most other 
points along the border of the hysteresis "ellipse" . From 
these inputs, we conclude that the phase boundary must 
run from the lower left to the upper right of the "ellipse" , 
intercepting the v = 4.5 line at 2.7T. Although the exact 
position of the phase boundary can not be derived from 
our data, a straight line as shown in Fig. ||d indicates its 
general position. Within the 4.4 < v < 4.6 band, the 
stripes align along the (110) direction to the left of the 
boundary. To the right of the boundary, the stripes align 
along the (110) direction. This analysis provides a phase 
diagram for the stripes around v = 9/2. Hystereses with 
much smaller areas are also observed around v — 11/2. 

It remains unclear why in previous studies samples 
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FIG. 3: Phase diagram of the orientation of the stripes at 
v = 9/2 in the B perp -B lp plane. Solid squares represent the 
(110) orientation of the stripes. Hollow triangles represent 
the (110) orientation. In the left panel, B lp points along the 
(110) direction. In the right panel, B !p points along the (110) 
direction. The inset shows the phase diagram at v — 11/2. 
Units are the same as in the u = 9/2 diagram. 

with densities inside the transition region of Fig. |3d nei- 
ther show the (110) orientation nor the hysteresis |l7|, p"8[ . 
It suggests the importance of the individual sample de- 
tails. Along this line, the thin AlAs layer at the interface 
may play a role in establishing the specific transition den- 
sity. 

The application of an in-plane magnetic field has be- 
come a useful tool in investigating the native symme- 
try breaking potential since it gives us additional control 
over the orientation of the stripes. In previous fixed den- 
sity samples, the stripes always aligned perpendicular to 
a sufficiently strong (^0.5T) in-plane field B lp [^2|, p3[ . 
This means that an in-plane field parallel to the B lp =0 
orientation of the stripes reorients the stripes to the per- 
pendicular alignment. 

We apply an in-plane magnetic field B lp to the 2DEG 
by rotating the HIGFET in situ in our dilution refrig- 
erator. Two cooldowns are required to align B lp along 
the (110) or the (110) direction respectively. The ro- 
tator assembly is able to reach a base temperature of 
70mK. Tilting angles were determined from a calibration 
of the rotator and from positions of well defined QHE 
and FQHE minima. The accuracy is better than 1°. By 
sweeping density at a fixed angle, i.e. fixed perpendicu- 
lar magnetic field, B perp , and fixed B ip , we determine the 
orientation of the stripes at v = 9/2. Fig. || shows the 
phase diagram in the B perp -B lp plane constructed from 
such data. The left panel represents data taken with B' p 
along the (110) direction. The right panel represents data 
taken with B lp along the (110) direction. A (110) orien- 
tation of the stripes is represented by a hollow triangle 
and a (110) orientation by a solid square. 



In order to assess the properties of the phase dia- 
gram, we first focus on the B vp —0 axis which simply re- 
flects the reorientation of the stripes discussed earlier: 
below BP er P=2.7T (equivalent to n= 2.9 x 10 u cm~ 2 , 
right axis), the stripes have the (110) orientation whereas 
above 2.7T, the (110) orientation prevails. As we in- 
crease B lp along either the (110) or the (110) direction, 
the (110) to (110) transition moves to higher densities. 
The density shift indicates a B lp induced energy gain 
of the (110) phase. It can be viewed as a sum of two 
terms: an isotropic term and an anisotropic term. The 
isotropic term is independent of the direction of B vp . The 
anisotropic term reflects the difference in energy gain be- 
tween a B lp parallel to the stripes and a B lp perpen- 
dicular to the stripes. Two conclusions can be drawn 
from the phase diagram. First, the roughly symmet- 
ric shape of the phase boundary indicates that in our 
sample, the dominant in-plane magnetic field effect is an 
isotropic energy gain of the (110) phase. The anisotropic 
term only becomes significant for \B lp \ > 0.4T. Second, 
in the high B lp regime, in the right panel, the transi- 
tion to the (110) phase occurs at lower densities than 
in the left panel. From this, we conclude that stripes 
oriented parallel to B w are energetically favored as com- 
pared to stripes oriented perpendicular to B lp . Both con- 
clusions are in contrast to previous fixed density sam- 
ples where the anisotropic term, which favors a per- 
pendicular alignment of the stripes with respect to B lp , 
dominates ]22], p3| , pi| , [2(|. Together with previous re- 
sults, we must conclude that the reaction of the stripes to 
an in-plane magnetic field is non-universal and sensitively 
depends on the details of individual samples. The unique 
density tunability of a HIGFET is essential to such a sys- 
tematic study. In fact, diverse behaviors are seen even 
within one specimen. As an inset to Fig. [| we show 
the phase diagram of the 11/2 state, constructed dur- 
ing the same cooldowns as the 9/2 state using the same 
method. The B ,p =0 transition also occurs at~2.7T. We 
can clearly see that the isotropic behavior observed at 
z/=9/2 is missing. The differences between the behav- 
iors of the 9/2 state and the 11/2 state strongly attest 
to the subtlety of the interactions involved. Examining 
other half- fillings, we note that the isotropic states of 
v = 5/2 and v — 7/2 become anisotropic in an in-plane 
field and the "hard" axis is always along the in-plane field 
direction. This is consistent with previous experimental 
findings and theoretical calculations [^2], pi[ . 

In summary, we have observed a density induced reori- 
entation of the stripe state at v = 9/2. As the density of 
the 2DEG is raised above 2.9 x 10 11 cm -2 , the stripes ro- 
tate from the (110) orientation to the (110) orientation. 
Our results demonstrate that the pinning direction of the 
stripes is not unique. 
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